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Abstract 

We investigate the presence of chaos in a system of two real scalar fields 
with discrete Z2 x Z2 symmetry. The potential that identify the system is 
defined with a real parameter r and presents distinct features for r > and for 
r < 0. For static field configurations, the system ntsupports two topological 
sectors for r > 0, and only one for r < 0. Under the assumption of spatially 
homogeneous fields, the system exhibts chaotic behavior almost everywhere in 
parameter space. In particular a more complex dynamics appears for r > 0; 
in this case chaos can decrease for increasing energy, a fact that is absent for 
r < 0. 
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I. INTRODUCTION 



Nonlinearity plays an important role in field theory. It is responsable for the presence 
of interactions and may enter the game allowing interesting situations, as is the case for 
instance when the system engenders spontaneous symmetry breaking. In this case the 
classical equations of motion may give rise to interesting field configurations such as kinks, 
vortices and monopoles, depending on the particular model in consideration. See for instance 
Refs. 0-0] for details. 

In the standard route to defect formation, one usually search for static field configurations 
that solve the equations of motion and present finite energy. In the simplest case of a 
single real scalar field, a system given in terms of a potential that presents Z2 symmetry 
may support topological solutions in 1 + 1 dimensions when the potential has at least two 
degenerate minima. In the case of two real scalar fields, the system is richer and may support 
distinct types of topological solutions, as we illustrate below. 

Nonlinearity also plays an important role on chaotic behavior of systems. For linear 
systems, the qualitative nature of the behavior does not change when one changes their 
parameters. However, for systems governed by nonlinear dynamics one may find examples 
where small change in a parameter can lead to dramatic changes in both the quantitative 
and qualitative behavior of the system. See for instance Refs. for details. 

There are distinct routes to chaos in field theory, and here we shall follow the point of 
view introduced in Ref. [0. In this case one investigates field theoretical models under the 
assumption of spatially homogeneous field configurations. We think of field configurations 
whose space variations are much smaller than the corresponding time variations, and so we 
treat the fields as depending only on time. This point of view has been explored in several 
different works, as for instance in Refs. [p|-[T0|] and in some works therein. We notice that 
these works deal with different models, describing Abelian and non Abelian gauge fields in 
the presence of spontaneous symmetry breaking. In these cases the symmetry to be broken 
is continuum, local, and the system supports vortices or monopoles. In the present work, 
however, we shall explore another system, simpler, that presents discrete symmetry and 
is described by a couple of real scalar fields. Recent investigations have shown that this 
system presents interesting properties and some soliton solutions have also been found, as 
we comment on below. 

The main motivation of the present work is to investigate the presence of chaos in a 
system that describes two real scalar fields engendering the discrete Z2 x Z2 symmetry. We 
shall do this by comparing the case where one supposes the field configurations to be spatially 
homogeneous with the more familiar situation which considers static configurations. The last 
case is appropriate for searching for soliton solutions, and we comment on that in the next 
Sec. |1[ We consider the case of spatially homogeneous fields in Sec [HI] where we search for 
chaotic behavior. As we also comment on, the subject of this paper is related to recent issues 



llj, which ask for instance whether the chaotic behavior found in flat spacetime persists 
during the cosmological expansion. In this paper we deal with a model which is defined in 
flat spacetime and may be seen as the results of an expanding FRW universe, valid under 
the approximation of very slow expansion rate. Furthermore, the present investigation lands 
very naturally to the context of hybrid inflation, where one requires real scalars such as the 
inflaton field and at least another scalar field, which couples to the inflaton field [F3| in a 
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way similar to the one we consider in this work. We end the paper by summarizing the 
results in Sec. ffV. 



II. GENERAL CONSIDERATIONS 

The system we study in this paper is a Z 2 x Z 2 model of two real fields t) and x{ x -> 
The potential that specifies the system is defined in terms of a single real parameter, r, 
which controls the number of minimum energy states of that potential. As we are going 
to show, when r changes sign the potential changes form, from the case with four minima 
(r > 0) to the case where only two minima are present (r < 0). For r > 0, spontaneous 
symmetry breaking appears in both the <fi and x directions in the (0, x) plane. For r < we 
have spontaneous symmetry breaking only in the 0-field direction, and this last situation is 
similar to the case of hybrid inflation where a first-order nonthermal phase transition after 
preheating seems to be present, as recently considered in Ref. [O. 

The field theory that we consider is described by the Lagrangian density 

C = \d a <\> «9> + \d aX d a x ~ x) (1) 

We are using standard notation, with h — c — 1, and V(<f>,x) is the potential, which one 
supposes to be given by 




mx) = 2(-) +o( — ) (2) 



where h = h(4>, x) is a smooth function of the two fields (f> and x- Here it obeys 

i . , . . r, . . n 
/' 



- h{(f>, x) = H<f>, X) = r Q</> 3 - <P) + <P X 2 (3) 



The parameter \x is real, with dimension of energy, and r is another real parameter, dimen- 
sionless. 



In 1 + 1 dimensions this system presents interesting soliton solutions [14], which has been 



used in applications in condensed matter [15[] and in field theory |I6}{18[|. The equations of 
motion for = <f>(x, t) and x = x{ x -> t) are given by 

<t>tt - 4>xx + V h H + h x h x<t> = (4) 

Xtt - Xxx + htj, hfa + h x h xx = (5) 

These are the equations we shall deal with in the following. For simplicity, however, we 
rewrite them in terms of dimensionless variables t = fit and x = fix to get 

4>e - 4>xx + hcf, h^> + h x h x<t> = (6) 

Xtt ~ Xxx + V + h x h xx = (7) 

In the case of static fields we have <fi = (f>(x) and x = x( x )- We change x — > y, for 
simplicity, and now the equations of motion become 
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= 2r 2 (0 2 -l)0 + 2(r + 2)0 X 2 (8) 
d 2 X 



dy 2 



2(x 2 -Ox + 2(r + 2)0 2 X (9) 



which are the equations we deal with when searching for soliton solutions. It is interesting 
to see that these equations of motion (|j) and @ can be solved by configurations that obey 
the pair of first-order differential equations 



__ rir _l)+^ (10) 

^ = 20x (11) 
dy 

Solutions that obey this pair of first-order equations are BPS solutions ||19|| . They are stable 
configurations that minimize the energy, as explicitly shown in |i"4|]15|| . See Ref. |18| for 



further comments on BPS solutions, and for showing explicitly that the system defined by 
the potential (|2|) is the bosonic portion of a supersymmetric theory. 

In the case of spatially homogeneous fields we have = (pit) and x — x(t)- We change 
t — > t, for simplicity, and here we get the equations of motion 



d 2 



dt 2 

fx 

dt 2 



2r 2 (l -0 2 )0-2(r + 2)0x 2 (12) 
2(r-x 2 )X"2(r + 2)0 2 x (13) 



These are the equations we have to deal with when searching for a chaotic behavior in the 
time evolution. It is interesting to see that the above equations have some analogies with 
the equations investigated in Ref. ||. However, in Ref. the system under consideration 
is the Abelian-Higgs model. This is the relativistic generalization of the Ginzburg-Landau 
theory of superconductivity, and it is well known that it presents vortex solutions [^(J . The 



reason for the similarity between eqs. (|12D and (|T3| ) and the equations investigated in |9fl 
seems to rely on the assumptions introduced in Ref. |J. 

We follow the aim of this paper, which is to investigate the chaotic behavior of eqs. ([12]) 
and ([HJ), together with the study of the static solutions of eqs. (|3|) and ([]). We start pre- 
senting some general considerations regarding the static solutions. Our system is identified 
by the potential V(4>, x) = L* 2 V(<p, x), where 

V(4>, X) = \r 2 (0 2 - l) 2 + r (0 2 - 1) x 2 + \ X 4 + 2 <f> 2 X 2 (14) 

The potential depends on the real parameter r, and then the behavior of the system is 
directly related to such parameter. The main characteristic of this potential, which certainly 
strongly affect both the static solutions and the chaotic time behavior, is its different shape 
for r positive and negative. In fact, when r > the potential contains four vacuum states 
in the (0, x) plane, while for r < there are only two vacuum states, and in principle we 
expect a different chaotic behavior for the dynamical trajectories that follows in accordance 
with eqs. (|l^) and (|T3|), as well as distinct features for the static solutions. It is important 
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to say that in the two particular cases when r = 1 and when r = —2 the two fields decouples 



17J and so we expect no chaotic behavior in these two cases. First we show that, from 



the point of view of the topological soliton solutions admitted by Eqs. (]8|) and (|9]), the case 
r > 0, r 7^ 1 is richer than the case of r < 0, r ^ —2. When r is positive, r ^ 1, the four 
vacuum states are given by 

vi = (-1, 0), v 2 = (1, 0) (15) 
U3 = (0,-VF), ^ 4 = (0,v^) (16) 

For r negative, r 7^ —2, we have the two vacuum states 

vr = (-1, 0), v 2 = (1, 0) (17) 

As it was already shown |I4| , |T5|| , the energy of the stable BPS solutions are obtained accord- 
ing to the values of h((p, \) in the vacuum states. 

For r positive, r 7^ 1 we name hi as the value of h at the vacuum state Vi, and so we have 

hi= (2/3) fir (18) 
h 2 = -(2/3) fir (19) 
^3 = ^4 = (20) 

This means that we have two type of BPS solutions, with energies Ei = \fi\ i = 1,2, where 

ei = ^r, (r>0, r^l) (21) 
e 2 = jjr, (r>0, r^l) (22) 

In the first case the solutions connect the vacuum states (±1, 0) and in the second case pair 
of vacuum states belonging to different axes. In the topological sector defined by the two 
vacuum states (—1,0) and (1,0) there are pairs of analytical solutions given by [0 

<t>{y) = ~ tanh[r(y - y)], x{v) = (23) 

and by 

0(l/) = -tanh[2( 2 /-y)] (24) 



Here y stands for the center of the topological solution. The first pair ([23] ) is valid for r > 
and the second one for r > 2. We see that the second pair of solutions (|24|) and (|25|) gets to 
the first pair ( p3|) in the limit r — >• 2. 

For r negative, r 7^ —2, there is just one type of BPS solution, with energy E\ = fx, 
where 

fx = ~|r|, (r<0, r^-2) (26) 
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The first pair of solutions given by Eq. (p3|) is also a pair of solutions in the case r < 0. 

It is clear that for r > 0, r ^ 1 the situation is richer than in the case r < 0,r^ —2, at 
least from the point of view of the above topological soliton solutions. In the first case the 
system admits two type of topological BPS solutions, while in the second case it has only 
one. We also expect a difference between the two r > and r < regimes to be present in 
the chaotic behavior of the eqs. fll2|) and ((Op. 

Another interesting feature of this model concerns the masses of the <ft and x fields. From 
the potential V(<p, x) we can introduce the matrix 



(27) 



where = d 2 V/d(j)d(f) and so forth. The superscript v stands for substituting the fields for 
their vacuum values once the derivatives are done. For the potential introduced in Eq. flUp 
we see that 

V 4>x = V^ = An 2 (r + 2)<Px (28) 

and so VL. = VjL vanishes for every vacuum state, despite the sign of r. This means that 
we can read the (square) mass of each field directly from the matrix (p?]) . We have at most 
two different mass values. They are 

m\ = 4 // 2 r 2 m 2 = 4 /j 2 (r > 0, <p axis) (29) 



x 

mi — 4 fj, 2 r m 2 =Ajj 2 r (r > 0, x axis) (30) 



and also 



m 



2 -4/xV m 2 =4/i 2 (r < 0,0 axis) (31) 



where we inform the axis of the vacuum state used to obtain the respective masses. We note 
that for r < the x field develops no symmetry breaking, and its mass does not depend on r. 
This is the way this specific Z 2 x Z 2 system behaves, and this behavior may have interesting 



connections with models of hybrid inflation [|T2yT3],^T]-g3| . We remark that there are other 



Z2 x Z2 systems, as for instance the ones investigated recently in Ref. [|I8" |, which are also 
defined with a single real parameter, but that present distinct behavior, with different sets 
of masses and vacuum states. They certainly lead to richer scenarios, and may give other 
informations on the role the parameter r plays in such models. 



III. CHAOTIC BEHAVIOR 

In this Section we investigate the chaotic behavior of the system introduced in Sec. [TI]. 
We focus our attention to the pair of Eqs. ( |12"D and (|i~3|). These equations can be seen as 
the equations of motion that follow from the Hamiltonian 

H = \vl + \vl + V{4>, X ) (32) 
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They describe the motion of a classical particle in the bidimensional potential V((f>, x)- We 
remark that this choice of Hamiltonian makes the energy dimensionless. The equations of 
motion can be written in first order form, in terms of the canonical variables (0, £></,) and 

(X,P X )> 

^ = wr n (33) 

dH 

p = -— = 2r 2 (l-0 2 )0 -2(r + 2)0 X 2 (34) 

X= dp~ X =PX 

Px = —Q- = 2 ( r - X 2 ) X - 2 (r + 2) 2 X (36) 

In the rest of the paper we study for which values of the parameter r and of the energy E 
the trajectories <f>(t), x(t) i Pcj>{t) i an d P x {t) show a regular or a chaotic behavior. 

The motion is always bounded, i.e. confined to a finite region whose size changes accord- 
ing to the energy E of the particle. In Fig. 1 we report the potential contour plot defined by 
x) = E\ the different values of the energy E are reported in the figure. Fig. la refers 
to r = —5 while Fig. lb and Fig. lc refer to r = 1 and r = 3, respectively. 

For r < we note the presence of two minima (±1, 0), corresponding to the two vacuum 
states previously discussed, and a saddle point at the origin. The potential vanishes at the 
two minima, while at the origin (0,0) it assumes the value r 2 /2 (12.5 for case r = —5 in 
Fig. la). Dynamical trajectories, e.g. solutions of Eqs. (j3~3~l)-([36|) are confined to one of 
the two minima regions if E < 12.5, while they can jump from one side to the other when 
E > 12.5. For r = —2 the two fields decouple and the numerical investigation shows no 
chaotic behavior. 

The situation for r > is qualitatively different because we have four minima V = at 
(±1, 0) and at (0, ±y/r), corresponding to the four vacuum states previously discussed, and 
a local maximum V = r 2 /2 at the origin (0, 0). There are also four saddle points with ^ 
and x 7^ 0; they are given by 



r 



2 



2(r + l)' 2(r + l) 

In the saddle points the potential gets the value 



r = (37) 



„2 



V = — r (38) 

2 r+1 y J 



The case r = 1 is a particular case because the potential presents symmetry |L7j (see 
Fig. lb); the two fields decouple and the numerical simulations show no chaotic behavior. 
When r ^ 1 trajectories can have different behaviors according to the energy E. We discuss 
the case r = 3 in Fig. lc: when E < 1.125 the four minima are unconnected regions and 
the motion remains confined around the minimum in which we start the trajectory. As soon 
as E > 1.125 (see Fig. lc) the trajectory crosses the saddle point [see Eq. (0)] between 
the minima and wonders from one minimum to another. The region around the relative 
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maximum at the origin (0,0) is still not allowed. When E > r 2 /2 = 4.5 the dynamical 
trajectories can cross the origin and the four mimina are all directly interconnected. 

In order to study the dynamical behavior of our system for different values of r and 
energy E, we integrate numerically Eqs. d33|) -(p6|) using a fourth order symplectic algorithm 
with a time step At = 0.001 p^fl . The time step has been determined in order to keep the 
error in energy conservation below AE/E = 10 -8 for any value of r and E; the results are 
stable respect to a further reduction of At. For any conservative system such as the one 
we are considering, the Hamiltonian is an integral of motion, and the energy conservation 
restricts trajectories to lie on a three-dimensional surface H(<fi,x,Pcf>,Px) = E in the four- 
dimensional phase space (<f>,XiP<f>iPx)- We can obtain a graphical information about the 
system by plotting the intersection of this three-dimensional surface with a plane. These 
kind of plot is called Poincare surface of section, and gives an indication about the dynamical 
behavior of the system. To define the surface of section in our case we follow a trajectory 
and we plot and x eacn time = 0. Regular regions will appear as a series of points (a 
mapping) which lie on a one dimensional curve (invariant KAM curve), while chaotic regions 
will appear as a scatter of points limited to a finite area due to energy conservation [Tj!| . 

In Fig. 2 and in Fig. 3 we show the surface of section in the plane (<j), x) for r = —5 and 
for r = 3, respectively. Each figure is obtained from 100 different trajectories followed for 
a time t = 500. The trajectories are generated from random uniformly-distributed initial 

conditions o ,JV o > Xo,P xo - 

For r = —5 we plot the two cases E = 10 and E = 20 (Fig. 2a and Fig. 2b respectively). 
When E = 10 the two minima are separated and the system exhibits only regular behavior 
(the surface of section shows only invariant KAM curves). When E = 20 the two minima 
are connected and the surface of section contains both regular and stochastic regions. The 
volume of the stochastic region increases with increasing energy, as it will be clear when we 
calculate Lyapunov exponents. 

For r = 3 the situation is richer. We report the surface of section for three different 
energies E = 0.5, E = 2, and E — 6. When E = 0.5 (Fig. 3a) the four mimina are not 
interconnected and we only have invariant KAM curves. When E = 2 (Fig. 3b) neighbour 
minima regions are connected. Regular and stochastic regions cohexist in the surface of 
section and are interwingled in a very complicated way. The scatter of points representing 
the chaotic region fill the area between regular curves, and then smaller regular islands are 
imbedded in the chaotic sea. The white regions correspond to regular region, and these would 
be filled by regular curves if we could increase even more the number of initial conditions 
considered in the figure. For energy E = 6 (Fig. 3c) the trajectory can cross the center 
and the resulting surface of section shows an increase of the chaotic area respect to the case 
E = 2. 

The cases r = 1 and r = —2 (which we do not report in figure) are two very particular 
cases. The two equations decouples, and indeed the Poincare surface of section show only 
KAM curves for every value of the energy E. 

The information we can get from Poincare surface of section is qualitative. A way to 
quantify the chaotic behavior of a system is by calculating its Lyapunov exponents. Chaos is 
defined in terms of the dyamical behavior of pairs of orbits which initially are close together 
in the phase space. The Lyapunov exponents are given by the rate of exponential divergence 
of close orbits and are defined from the long term evolution of an initial infinitesimal volume 
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f2 in the phase space from the following formula 

A; = lim XAt) = lim - In (39) 

t-*oc iV ' t^oo t di(0) V ' 

where di(t) are the lenght of the principal axes of the ellipsoid Q at time t. There is one 
Lyapunov exponent for each dimension of the phase space and the flow is chaotic if 
at least one Lyapunov exponent is positive. For a two-dimensional Hamiltonian system we 
expect two positive Ai, A2 and two negative A3, A4 Lyapunov exponents with the following 
symmetry A3 = — A2 and A4 = — Ai, because of the symplectic structure of Hamiltonian 
systems. 

To compute Lyapunov exponents for different values of r and energy E we use the 



standard method developped in Ref . |25[ . The average number of time steps to use in order 
to have a good convergence of the results is about 10 7 — 10 8 , which means t = 10 4 — 10 5 . The 
symmetries A 3 = — A 2 and A 4 = — Ai are always perfectly verified in the numerical results; 
this is a good check for the integration codes. Fig. 4 shows the numerical calculation of the 
two largest Lyapunov exponents for case r = 3 and E = 5.5; we report the quantity Xi(t) for 
i = 1, 2 as a function of time. The top panel in Fig. 4 shows two different trajectories both 
starting from initial conditions in the chaotic region. After a transient time, X\(t) converges 
for both trajectories to a constant value of about 0.3, while X 2 (t) converges to 0.01. To 
obtain A, we evaluate numerically the limit in formula ( |39"D considering a time average of 
Aj(t) over the time interval t = 2 ■ 10 4 — 3 ■ 10 4 . The bottom panel of Fig. 4 shows the 
behavior of Ai(£) and A 2 (t) for a different initial condition, this time chosen to be in the 
regular region of the phase space; Xi(t) and A 2 (t) do not reach a finite asymptotic value but 
tend to zero for large t, showing that the Lyapunov exponents are equal to zero. 

For each value of r and energy we consider N = 100 different trajectories with random 
initial conditions 0o, p</> , Xo, Vxq- We compute Aj for each trajectory and we count the number 
N c of chaotic trajectories (i.e. with positive Ai, as in Fig. 4a). In the following we report the 
ratio R = N c /N as an indicator of the measure of volume of phase space which is chaotic, 
and the value of Ai averaged over the N c chaotic trajectories. We do not report A 2 which is, 
for each value of r and energy, one order of magnitude smaller than Ai. 

We consider separately the two cases r < and r > for which we expect to find a 



qualitative difference in the chaotic properties of Eqs. (|i2|) and (|13|). We discussed above 
the topological soliton solutions of Eqs. (^) and (0) for these two cases and we found that case 
r > is richer than case r < 0. In fact, for r > 0, r 7^ 1 the system supports two topological 
sectors, admitting two different type of topological BPS solutions; for r < 0, r 7^ —2 there is 
just one topological sector. This difference is due to the presence of four or two minima in 
the potential V(<f>, x), for r > or r < respectively. We expect the difference in the shape 
of the potential to have important implications also on the chaotic behavior of our system. 
In Fig. 5 and Fig. 6 we plot the ratio R = N c /N and Ai, as function of energy for the two 
cases r < 0,r ^ —2 and r > 0,r ^ 1, respectively. The system exhibits an order to chaos 
transition as a function of energy for any value of r, but for r = —2 and r = 1. These two 
cases are particular cases because the two field decouple and the system degenerates to two 
systems of a single field each one; R and Ai are identically equal to zero in the whole energy 
range and are not reported in figures. 

We show the typical behavior for r negative in Fig. 5. We report R and Ai vs. energy 
for r = —5; other values of r have the same qualitative behavior. The onset of chaos is for 
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energy equal to r 2 /2 (12.5 in figure), when both minima can be visited by a single trajectory 
and R and Ai start to be different from zero. For larger values of energy both R and Ax 
increase with the energy. 

The situation is different for r positive. In Fig. 6 we plot R and Ai vs. E for two different 
values of r, r = 1.5 (dashed line) and r = 3 (dashed-dotted line). The behavior for r = 1.5 
and r = 3 is qualitatively similar, although shifted in energy. Chaos starts when the energy 
overcomes the barrier between two minima [which is given in accordance with Eq. (|58D], 
and it keep increasing as function of the energy until a new back bending appears, both in 
R and in Ai. This behavior is to be connected to a stabilization effect occurring when the 
trajectory can cross the origin since both R and Ai have local minima at E = r 2 /2. Although 
the chaotic behavior for r = 1.5 and for r = 3 are qualitatively similar, the absolute value of 
R and Ai is as larger as closer r gets to unity. This result can be of interest to application 
since the absolute value of r also controls the masses of the two fields and may be used in 
models of hybrid inflation. 

The behavior at r > is therefore different from the behavior at r < 0. The most 
important qualitative difference between these two cases is that for r > we can find 
localized regions in energy where suppression of chaotic behavior appears for increasing 
energies. Such behavior is not present for r negative, and this is directly related to the 
reduction of the number of minima in this last case. Similarly, we have shown that for static 
solutions the number of topological sectors changes from two to one when r changes from 
positive to negative values, respectively. 



IV. COMMENTS AND CONCLUSIONS 

In this work we have investigated the presence of chaotic behavior in a system of two 
real scalar fields that engenders the Z 2 X Z 2 symmetry. The system is defined by a potential 
that is controlled by a single real parameter, r. This parameter can be positive or negative, 
and each case leads to different behavior. For static field configurations for instance, the 
system suports two dintinct topological sectors for r positive, r ^ 1, and only one when r is 
negative, r ^ -2. 

Within the context of chaotic behavior, we have shown that chaos is present almost 
everywhere in parameter space, with distinct qualitative behavior for r > and for r < 0. 
The system is richer for r positive, and this is directly related to the presence of the four 
minima when r > 0. For r > we can find localized regions in energy where suppression of 
chaotic behavior appears for increasing energy values, a fact that is absent for r negative. 
On the other hand, in the case of static solution there are two topological sectors for r > 0, 
in contraposition with the single sector that the system supports for r < 0. 

There are other models that engender the Z 2 x Z 2 symmetry, and that are also governed by 



a similar real parameter r. Some examples appear in Ref. flSfl , and they can also be studied 
within the context of chaotic behavior, under the assumption of spatially homogeneous field 
configurations. Such investigations can introduce further light toward a better understanding 
of the connection between chaotic behaviors and the parameter r that specifies the physical 
properties of the system. 

As we have already commented on, the present investigations may be of interest to infla- 
tionary cosmology, in connection with issues raised in the recent works pT ^T3] , pT| -p3]] . The 
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present results seem to be valid in an expanding FRW universe when the rate of expansion 
is very small. The chaotic behavior that we have found may have been present in early times 
and may have played some role in the cosmic evolution. A study in which the coupling to 
gravity is fully considered, and other related issues are presently under cosideration. 
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FIGURE CAPTIONS 

Fig. 1. Potential contours plot for different values of E. a), b), and c) refer to r = —5, 1 
and 3, respectively. 

Fig. 2. Poincare surface of section for r = —5 in (0, \) plane. E = 10 in a) and E = 20 in 
b). 

Fig. 3. Poincare surface of section for r = 3 in ((f), x) plane. E = 0.5 in a), E — 2 in b), 
and E = 6 in c). 

Fig. 4. Lyapunov Exponents Xi(t) and A 2 (i) vs. time for r = 3 and E = 5.5. Two chaotic 
trajectories are plotted in the top panel, while a regular trajectory is plotted in the bottom 
panel. 

Fig. 5. Fraction R of the phase space which is chaotic and largest Lyapunov exponent Ai 
vs. E for r = —5 

Fig. 6. Same as in Fig.5 for two positive values of r, namely r = 1.5 (dashed) and r = 3 
(dotted-dashed). 
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